Unveiling the Landau Levels Structure of Graphene Nanoribbons 
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Magnetotransport measurements are performed in ultraclean (litiiograpiiically patterned) 
graphene nanoribbons down to 70 nm. At high magnetic fields, a fragmentation of the electronic 
spectrum into a Landau levels pattern with unusual features is unveiled. The singular Landau spec- 
trum reveals large magneto-oscillations of the Fermi energy and valley degeneracy lifting. Quantum 
simulations suggest some disorder threshold at the origin of mixing between opposite chiral magnetic 
edge states and disappearance of quantum Hall effect. 
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PACS numbers: 72.80.Vp,75.47.-m,73.22.Pr 



Introduction.- To benefit from the unusual trans- 
port properties of graphene (massless Dirac fermions 
physics) [ll, |2| , and engineer novel building blocks for fu- 
ture carbon-based nanoelectronics the fabrication of 
clean materials has become a central issue. Of great con- 
cern is the design of graphene nanoribbons (GNRs) with 
lateral sizes in the order of a few to tens of nanome- 
ters, which allow some gap engineering The trans- 
verse confinement leads to ID electronic sub-bands pro- 
files whose details depend on the width and the edge 
geometry of the ribbons 

In presence of a strong perpendicular magnetic 
field, anomalous quantum Hall effect develops in two- 
dimensional graphene with unique properties widely dis- 
cussed in the literature The electronic spectrum is 
theoretically predicted to evolve into magneto-electronic 
sub-bands resulting from a competition between mag- 
netic and electronic confinement jil, Hj. This is partly un- 
veiled by anomalous Shubnikov-de Haas oscillations when 
the Landau diameter becomes larger than the ribbon's 
width However, it is puzzling to note the lack of ex- 
perimental evidence of Hall quantization in nanoribbons 
narrower than 200 nm Recent magnetotransport 

experiments in chemically derived narrow GNRs reported 
some signatures of chiral magnetic edge states revealed by 
a large positive magnetoconductance [l^l ■ Positive mag- 
netoconductance is also observed in lithographic nanorib- 
bons [ll| . but in all these experiments, the conduc- 
tance is far from being quantized, and transport remains 
strongly diffusive, thus jeopardizing a convincing obser- 
vation of the underlying Landau levels. Several sources 
of disorder are suspected to crosslink chiral edge currents 
for the narrowest ribbons, thus preventing quantum Hall 
effect from developing Additionally, for clean GNRs, 
recent theoretical calculations suggest a strong impact 
of electron-electron interactions on the band structure 



leading to the suppression of the conductance quantiza- 
tion [l3| . A detailed characterization of Landau levels in 
GNRs thus remains to be accomplished. 

In this Letter, we report two-terminal quantum Hall 
resistance measurements on GNR devices. The oscilla- 
tory behavior of the magnetoresistance (MR) evidences 
a clear signature of magneto-electronic subbands. The 
onset of a singular graphene-like Landau spectrum re- 
veals a large pinning of the Fermi energy along with a 
valley degeneracy lifting. To rationalize those features, 
we simulate the spatial extension of the corresponding 
magnetic edge states and their distribution in presence 
of disorder. 
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FIG. 1: (color online) Experimental G(Vg) curves measured 
at 4.2K on two GNR devices of width 100 and 70nm (sample 
A and B). Inset : the AFM image of the devices. 



Graphene devices are made by mechanical exfoliation 
of graphite onto 7T,++Si/SiO2(300nm) substrates. Elec- 
trodes are defined by e-beam lithography followed by 



2 



thermal evaporation of Ti(lnm)/Pd(10nm)/Au(40nm). 
Next, the graphene flakes are patterned into ribbons by 
oxygen plasma etching, using PMMA as an etching mask. 
A set of connected GNRs are prepared with width {W) 
ranging from 60 to lOOnm and length (L) from 350 to 
800nm. In the following, we present extended results 
on two devices (inset Fig.l), sample A (L = 350nm, 
W = lOOnm) and B (L = 750nm, W = 70nm), having 
the hallmarks of the overall samples. 
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FIG. 3: (color online) (a) Experimental magnetoconductance 
of sample A versus the filling factor deduced from the MR 
curves at selected Vg, from -40V (top) to OV (bottom), by 
step of 5V. The dashed line is the simulated G{i').{h) Magne- 
toconductance of sample B versus the filling factor at different 
from -50V (top) to OV (bottom), by step of 5V. Inset: the 
singular MR curves exhibiting a double resistance peak when 
Ef crosses the A'^ =2 Landau level. 



FIG. 2: (color online) Two-probe perpendicular magnetore- 
sistance measured at 4.2K on sample A, for selected Vg from 
-lOV (cyan) to -40V (light magenta), by step of 5V. The inset 
shows anomalous 1/B Shubnikov-de Haas oscillations with 
circle marks indicating the Landau index A'^ as a function of 
the 1/Bi locations deduced from the MR curve at -40V. The 
red crosses are the N,{l/Bi) simulated data from the band 
structure (Fig. [4k, inset). The blue curve is the calculated 
N{l/Bi) from 



After standard thermal annealing treatment in vacuum 
to desorb contaminants, the measured G(V^) at 4.2K 
(Fig. [T]) exhibits a minimum at the charge neutrality 
point corresponding to a back-gate potential Vg = Vcnp , 
close to zero, 2.5V and —0.5V for samples A and B, re- 
spectively, thus pinpointing a negligible residual doping. 
From a numerical calculation of the electrostatic cou- 
pling between the ribbons and the back-gate (available 
at http://www.fastfieldsolvers.comm), the carrier den- 
sity is estimated as n(m~^) 



1.5 X IQi^ X (Vg(V) - 
Vcnp) while the Fermi energy scales as Epipi^) w 
40 X \J{Vg^{y) — Vcnp)- Assuming a negligible contribu- 
tion of quantum interferences, we estimate the electronic 
mean free path, lm{Vg), from the G{Vg) curves At 
Vg = -20V, we infer sa 50 - 120nm (80-120nm) for 
sample A(B). In both cases, Z,„ ^ W and L/l,n ~ 2 — 7, 
this indicates that the transport regime is close to a quasi- 
ballistic regime. The good quality of the samples is also 
confirmed by the estimated large field effect mobility at 
4.2K /XA(B) « 1200(3500)cm7(Vs), and by the Fabry- 
Perot conductance modulations observed at 2K [l3 |. 



Under a perpendicular magnetic field, the two-terminal 
resistance shows a non-trivial sample-shape dependent 
profile. It reveals fingerprints of both the quantized Hall 
resistance for incompressible charge carriers densities (i.e. 
for filling factors vijieh/ eB) = 4(n + l/2)) and the longi- 
tudinal resistance for intermediate filling factors [l5|. Ac- 
cordingly, for elongated ribbons, the resistance peaks pin- 
point the depopulation onset of the highest occupied Lan- 
dau subband of index N . Figure [5] shows the two-probe 
MR up to 55T at various Vg for sample A. An oscillatory 
behavior of the resistance along with quantized minima 
(/i/6e^) and resistance plateaus (/i/2e^) are clearly visi- 
ble. The MR plot as a function of the inverse magnetic 
field shows the expected 1/B Shubnikov-de Haas oscilla- 
tions of the resistance peaks for a 2D electron gas jl^ . 

However, a strong departure from the 1/B periodic- 
ity is observed for larger TV (inset. Fig. [2] circle marks). 
This is a convincing signature of the electronic confine- 
ment that starts to overcome the magnetic one when 
the cyclotron radius becomes larger than W/2. The lin- 
earity and its deviation above N ^ 9 are well repro- 
duced by the calculation of the number of occupied sub- 
bands N = N{l/Bi) in the frame of semiclassical Bohr- 
Sommerfeld quantization rule with an hard-wall confine- 
ment (blue line) [l4| . 

Figure [3] shows the conductance profiles of the two 
GNRs as a function of the filling factor. The curves are 
deduced from the MR experiments, i.e. a constant charge 
density along with a pulse field sweep, meaning an in- 
crease of the Landau energy broadening versus v. For 
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FIG. 4: (color online) (a) and (b) Band structure of a 70nm- 
wide aGNR around the CNP at and SOT, respectively, (c- 
e) Simulated magneto-electronic subbands versus magnetic 
field for the lOOnm-wide aGNR (c) and 70nm-wide aGNR 
(d,e) directly compared to quantum oscillations at selected 
Ef- Black lines hold for the edge subbands at zero-fc. The 
dashed black Unas correspond to the valley degeneracy lift- 
ing at nonzero k. In red are plotted the MR curves in arbi- 
trary unit and in blue, the related Fermi energy. Dark yellow 
vertical lines indicate the resistance peaks shifted to higher 
magnetic field due to the Fermi energy pinning at low N . 



sample A (Fig. 2a), a clear 2Go Hall conductance plateau 
is observed at = 2, providing direct evidence of a single 
layer graphene. For larger incompressible charge densi- 
ties, at 1/ = 6, 10, only maxima of conductance develop 
instead of the expected (6, 10, ...)Go plateaus. This over- 
all behavior is consistently described in the context of 
the conformal invariance of the conductance, where the 
distortion of the plateau is driven by the device aspect 
ratio ^ = LjW . Following the theory developed in [15|, 
we reasonably well simulate our data for t,fit = 4.1 (to be 
compared to ^ex^ = 3.5) and a broadening of /\v sa 0.84 
and 1.05 around v = 2 and 6, respectively (Figl3^, dashed 
line). The larger broadening at higher filling factor is 
therefore responsible for the shrinkage of the 6Go quan- 
tized state. 

Despite such agreement, several intriguing experimen- 
tal features demand further considerations : (i) The re- 



sistance peaks preceding the hj^e^ plateau on sample A 
(i.e. corresponding to the crossing of the iV = 3 Landau 
level) appear drastically enlarged for the MR curves at 
larger (Fig[2]). (ii) Surprisingly, the Hall conductance 
of sample B, which is 30% narrower, does not present 
a well-defined 2Go plateau even though the 6Go quanti- 
zation at = 6 is preserved (Figl3l3). The suppression 
of the 2Go conductance plateau goes along with a well- 
marked splitting and broadening of the resistance peak 
that develops before the expected plateau (Figl5]D-inset). 
(Hi) Both samples exhibit a gradual suppression of the 
6Go conductance at j/=6 as Vg decreases (Figl3^ and b) . 

For a deeper understanding of the MR curves and the 
related Landau levels pattern, we now consider the mag- 
netic field dependent band structure of two armchair rib- 
bons (aGNRs) of width W k. 100 and W « 70nm. Due 
to the rather large width of the ribbons, the following dis- 
cussion does not depend on the exact number of dimer 
lines that compose the ribbons. Figures |4l[ a, b) show the 
band structure of the narrowest aGNR around the CNP 
at i3=0 and 50T, respectively. Note that the valley de- 
generacy of the two-dimensional graphene is lifted as a 
consequence of the ribbon boundary conditions At 
B=50T, the bulk Landau levels are degenerated while the 
edge states (at large wavenumber /c) , split into couples of 
nondegenerate bands, among which one has an increasing 
energy while the other first decreases and then rises when 
approaching the edges of the Brillouin zone. The black 
(solid and dashed) lines of Fig. Hljc-e) indicate the mini- 
mum energy of each band as a function of the magnetic 
field. Above lOT, they start to scale as Vs, thus con- 
verging to the Landau levels scaling of the 2-D graphene. 
Interestingly, the degeneracy breaking of edge states is 
enhanced at high fields. To relate the transport oscil- 
lations with the underlying band structure, we compare 
the intersection of the Fermi energy and the magneto- 
subbands spectrum with the locations of the maxima of 
resistance (red curve). Assuming a constant Ep (blue 
dashed line), a good agreement between the locations of 
the resistance peaks and the subband depopulation for 
high quantum numbers is observed at low fields (inset 
of Fig. and corresponding red crosses in the inset of 
Fig. [2]). However, a noticeable mismatch gradually in- 
creases below A^=3 and reaches several Tesla for = 2 
in Fig. IDJd.e). marked by vertical dashed lines [l^ . 

To shed light on this issue, we observe that, as in 
conventional 2D gas, the Fermi energy is not constant 
but varies to accommodate the carriers into the avail- 
able magnetic field dependent subbands (lj|. This de- 
termines large oscillations of the Fermi energy at low N 
(see the blue lines in Fig. U]) , with an evident pinning at 
the Landau levels at high magnetic fields. The Fermi en- 
ergy oscillations provide a clear explanation of the shape 
of the MR curves in the high field regime. Indeed, the 
inflexion point of the Ef(B) curves in between two suc- 
cessive Landau levels well matches with the minima of 
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resistance (Fig. |3)Ic,d); marked by arrows). Besides, the 
widening of the resistance peaks at larger Vg and their 
shift to high fields are a direct consequence of a stronger 
Fermi energy pinning onto the lowest index Landau lev- 
els • In this framework, the absence for sample B of a 
well-defined 2Go plateau accompanying the splitting and 
the widening of the resistance peak also finds a natural 
explanation (Fig. and [3|)- The second maximum 
of resistance coincides with the changeover oi Ep from 
the N ~ 2 Landau level (solid black line) to the sec- 
ond energy minimum at nonzero k (dashed black line). 
Note that such a signature of the valley degeneracy lifting 
on the MR curves requires a large enough energy split- 
ting between the two sub-levels along with a pinning of 
the Fermi energy on the two states. These conditions 
are borne out on sample B, which has a higher mobility 
and allows for a sharper valley degeneracy lifting with an 
higher density of states at the subband edges. 
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starts to degrade when the corresponding edge channels 
in the pristine ribbon are still well separated by 40 and 
15nm for the sample A and B, respectively. This is con- 
sistent with the higher degree of disorder of sample A, 
which is able to crosslinks the chiral currents more effec- 
tively. This behavior is very sample-dependent, since a 
specific disorder configuration might make the chiral cur- 
rents come significantly closer, thus facilitating backscat- 
tering. Additionally, the lower Vg the more pronounced 
the phenomenon is, owing to a reduced screening of the 
impurity potential. As a matter of illustration. Fig. 
shows the leakage of edge currents (at = 6) provoked 
by a single charged impurity for a 70nm-wide graphene 
ribbon. 

Conclusion.- Unusual features of Landau levels have 
been explored in clean graphene nanoribbons using high 
magnetic field transport measurements. The observed in- 
triguing transport oscillations have been related to Fermi 
level pinning and valley degeneracy lifting. By means of 
quantum simulations, the spatial extension and robust- 
ness to disorder of chiral edge currents have been dis- 
cussed. 
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FIG. 5: (color online) (a) The measured suppression of the 
quantized conductance at = 6 for sample A (red) and B 
(black) versus the simulated current-free bulk width AZ(nm). 
(b) Spatial distribution of the background spectral currents 
through a 70nm-wide aGNR at = 6 (Ep = 119meV and 
B = 12T), in presence of a Gaussian potential with maximum 
strength V = 1 eV, mimicking a single charged impurity. 



We finally comment on the gradual suppression of the 
6Gq conductance at = 6 when decreasing the back- 
gate voltage (Fig. The deterioration of the conduc- 
tance quantization at an incompressible filling factor can 
be related to the presence of disorder that crosslinks the 

14| , we sim- 



chiral currents flowing at opposite edges. In 
ulate the spatial profile, along the transverse section of a 
pristine 70nm-wide aGNR, of the two counter propagat- 
ing edge channels at ;/ = 6 as a function of the magnetic 
field. As the magnetic field increases, the Landau states 
extend more and more in the bulk and the edge states 
are progressively pushed toward the edges. We extract 
the distance between the edge channels (A/) for a given 
magnetic field. Figure [S^ shows the measured conduc- 
tance at = 6 as a function of Al. The conductance 
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